One- Way Hysteresis in the Depinning Transition 
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We identify a mechanism for a type of hysteresis which we predict to occur in a variety of depinning 
transitions. We show that the phenomenon of one-way hysteresis is generic to stress-overshoot 
models of the depinning transition, and we show how to calculate the size of the hysteresis gap 
exactly in a large class of models. Unlike thermodynamic phase transitions, the dynamic phase 
transition is still continuous despite the presence of hysteresis because the terms which produce 
the hysteresis are renormalization group irrelevant. We discuss the experimental and numerical 
signatures of one-way hysteresis, which includes microscopic nucleation. 



Many dynamical phase transitions display hysteresis in 
the transition from the static to the moving phase. The 
most familiar of these is solid-on-solid friction- the value 
of the applied force which causes a solid to start sliding is 
greater than the force at which a moving solid will stop. 
Another well studied example is the propagation of cracks 
in stressed solids. The crack moves forward at a nonzero 
velocity when the force on the solid exceeds a critical 
value 0] , and as the force is reduced, it is likely that the 
crack stops moving at a second critical force lower than 
the first. However, the crack front posesses a well-defined 
roughness exponent in a variety of experiments Q , which 
suggests that the depinning transition is continuous. 

Hysteresis also occurs in other dynamical depinning 
transitions transitions: the current-flow transition in 
charge-density wave solids 0, the depinning of vortex- 
lines in type II superconductors Q , the onset of motion 
in hctereogeneous faults 0, and the wetting of a rough 
surface by helium-4 6J. Previously it had always been 
assumed that the hysteresis implies that the transition 
is discontinuous, by analogy with the world of thermo- 
dynamic phase transitions where hysteresis only occurs 
in the first-order ones. We show that in noncquilibrium 
phase transitions this is not so, and that a particular one- 
way hysteresis of the type observed naturally occurs in 
continuous phase transitions. This suggests that, among 
the systems we observe to have one-way hysteresis, a 
number of them can still be described by a statistical 
field theory at their critical point. We can show this be- 
cause the perturbation that generates one-way hysteresis 
is renormalization group irrelevant, so the models with 
hysteresis remain in the same universality class with the 
same critical exponents as in the absence of the pertur- 
bation. The universal region is in a hard-to-reach corner 
of parameter space, however. To get there, we must first 
increase the applied force to start the system moving, 
and then decrease it until the velocity is small. If a large 
fluctuation causes the system to stop, we must repeat the 
procedure. This makes both numerical and experimental 
measurements of transition quantities arduous. 

The source of the hysteresis in our model is the nonad- 
ditive component of the transient stresses in the system, 



and we can subtract it out by a simple procedure. The 
subtracted models are equivalent to the unsubtracted 
ones, not only in the universal region but everywhere, 
except that the one-way hysteresis is removed. The sub- 
traction procedure is the main calculation tool, and it 
allows us to demonstrate that the hysteresis gap exists, 
and to calculate the size of the gap in our models exactly. 
Performing the subtraction makes numerical analysis of 
the transition much simpler, and the subtracted models 
allow us to investigate the phase diagram of the contin- 
uum limit more easily. 

The models we analyze are minor generalizations of 
the stress-overshoot model introduced by Schwarz and 
Fisher 0, Earlier, Fisher et al. had suggested 
that models of this sort might possess a nonvanishing 
hysteresis gap. Surprisingly, for the particular stress- 
overshoot model in this is not so — the hysteresis 
vanishes in the infinite system limit. But the hysteresis 
is abnormally persistent, and this indicates that a small 
modification should produce true hysteretic behavior. 

Our models have a finite hysteresis gap in the infinite 
system limit, proving Fisher's conjecture. Furthermore, 
the one-way hysteresis occurs in the generic case — ear- 
lier models had been accidentally fine-tuned so that the 
hysteresis was a finite-size effect. But the existence of 
the gap is paradoxical- the term that produces it is irrel- 
evant in the renormalization group sense, and should not 
affect the critical behavior. We resolve the contradiction 
by showing precisely how an irrelevant perturbation can 
introduce one-way hysteresis without otherwise altering 
the critical behavior [10( . The close relation between the 
stress-overshoot models and our models provides an ex- 
planation for the persistence of the finite-size hysteresis 
in y|. 

The model describes a surface whose position /i(x) in- 
creases by jumps. The position of the surface describes 
the contact of two solids, or the position of a crack front, 
or the phase in a charge-density wave. There are pinning 
forces which hold the surface back, and a global driv- 
ing force which pushes the surface forward. When the 
total force on a site is greater than the pinning force, 
the site jumps. The sites communicate via an elastic 
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force, which is linear, as if every site was connected to 
its nearest neighbor by springs. In addition, there are 
transient velocity-dependent forces on neighboring sites 
which vanish after a short period of time. These stresses 
are analogous to inertial terms in highly dissipative equa- 
tions, and like those terms they are generically irrelevant. 
These are the stress overshoots which are responsible for 
the hysteresis. 

The stress overshoots model a variety of very different 
physical effects. In a propagating crack front, the stress 
overshoots model the effect of the acoustic wave which 
propagates in the solid whenever a segment of the crack 
moves forward. This acoustic wave pushes on nearby 
segments, possibly causing them to depin. In models of 
solid-on-solid friction, a stress overshoot can model the 
transient component of the adhesive force between the 
solids, which rises to a maximum value over a period of 
time. 

To study these effects, on every site x of a d- 
dimensional lattice define: 

h(x) - the position (height) of a d-dimensional surface, 
initially zero. 

w(x) - the velocity, either zero or one. 

/(x) - a pinning force uniformly random in [0, P]. 

There are three global parameters: 
P- the maximum strength of the pinning. P > 2d to 
avoid lattice artifacts. 
F- the global applied force 
M— the magnitude of the stress overshoots. 

To update the model: 

1. Calculate the elastic- stress E as the sum of /i(y) — /i(x) 
over all nearest neighbors y. 

2. Calculate the overshoot- stress S as the sum of v(y) 
over all the nearest neighbors y including x itself in the 
sum 

3. At every site xwhere E(x) + MS(x) + F > /(x), set 
v to one, elsewhere set it to zero. 

4. Add v to h at every site. Wherever v is nonzero, the 
site is said to have hopped. 

5. At every x that hopped, generate a new pinning force 
/in[0,P]. 

This model is a discrete form of a standard model for 
the depinning transition. The parameter we vary is the 
applied force F, and the order-parameter is the average 
velocity v. When no overshoots are present (M = 0), it 
has a continuous transition at the critical force F c . This 
transition has been extensively studied, and the critical 
exponents have been calculated to first- 0,0,0 and 
second- |l5j order in epsilon. 

The key to understanding the transition in the M = 
model is the no-passing rule of Middleton 0] . To state 
the no-passing rule, we must imagine that the random 
forces / have been pregenerated at the beginning of the 



simulation, and they are discovered as the surface moves 
forward. The pregenerated force may be thought of as 
a random function of the height h(x). The no-passing 
rule states that when two surfaces fti(x) and /i2(x) move 
in the same pregenerated pinning, the surface which is 
initially behind will at no time overtake the surface that 
is initially ahead. 

The proof is by induction in the time. Assuming that 
h\ is ahead of hi at time t, in order to determine whether 
h\ is still ahead at time t + 1 we only need to examine 
the points where hi and hi coincide. At these points, the 
pinning force on h\ and h^ are equal because the pinning 
is a function of the height. The elastic force on h\ is 
greater than or equal to the elastic force on hi because 
the elastic force is an increasing function of the position of 
the neighbors. Therefore, wherever hi jumps, h\ jumps 
too, and hi is not ahead at time t + 1. 

Only assuming that the steady-state behavior of the 
surface is motion with a constant velocity, the no-passing 
rule allows us to conclude that the velocity must be a 
single- valued function of the applied force F, since at any 
force where two velocities may coexist, a surface moving 
with the faster velocity will overtake a surface moving 
with the slower velocity, even if it is initially behind. 
Since it is intuitively clear that the velocity is an increas- 
ing function of the applied force, the no-passing rule also 
makes it all-but-certain that the transition in this model 
is continuous, since at any jump-discontinuity the limits 
v(F+) and v(F-) are different, and we may start two 
surfaces at two infinitesimally different forces and let the 
faster overtake the slower 0|. While these arguments 
as stated are not mathematically rigorous, they are sig- 
nificant, because this transition is one of the few whose 
order may be established theoretically. 

From the proof, it is clear that no-passing is a generic 
feature of any elastic force law for which the elastic force 
increases when the height of the neighbors increases. 
This implies that adjusting the elasticity law in any rea- 
sonable way cannot produce hysteresis — a new physi- 
cal mechanism is needed. Because the stress overshoots 
are transient, they can circumvent no-passing. A mov- 
ing surface can move past a stationary one because the 
stationary one feels no overshoots. 

Now we will demonstrate that the addition of stress 
overshoots of the type we have considered gives rise to 
hysteresis. To make the proof transparent, we will intro- 
duce three simple models, one of which has no hysteresis, 
and the other two which have a one-way hysteresis gap of 
exactly M units of force. We will show that they are all 
related by a subtraction procedure that removes M units 
of one-way hysteresis. Applying the subtraction proce- 
dure to our model we will arrive at a subtracted model 
which possesses a hysteresis gap which is exactly M units 
of force smaller. Irrelevance of the perturbation implies 
that there is no hysteresis at all in the subtracted model, 
from which we conclude that our model has exactly M 
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units of hysteresis. 

The first model we introduce has a stress overshoot 
which makes it trivial to calculate the hysteresis gap. The 
model is defined with an overshoot stress S(x) which is 
a global OR of all the Boolean velocities, that is, it is 
equal to 1 if there is a single nonzero velocity anywhere 
in the lattice, and it is equal to zero only if the surface is 
motionless. This is a global perturbation of the model, 
and so we have no right to expect that the behavior is 
unchanged. But the behavior of this model at any M is 
obvious — the overshoots simply shift the critical force 
by M units downward when the interface is moving, and 
they do nothing at all when the interface is static. The 
effect of the global overshoot is trivial- the model will 
start moving at a critical force F c at a finite velocity, 
and it will stop moving at a critical force F c — M where 
it will display identical behavior (site by site and hop 
by hop) to the M — model. The transition remains 
continuous and the critical exponents are unchanged. 

The second model we introduce is local, and is superfi- 
cially different from the first. In this model, which we call 
the nonadditive model, the overshoot stress is a local OR 
of the neighboring sites' velocities- S(x) = 1 only if v(x) 
is nonzero or v(y) is nonzero for one of the nearest neigh- 
bors y of x. But difference is superficial- the nonadditive 
model is (site by site and hop by hop) equivalent in the 
steady state to the global model. In particular, it has a 
continuous phase transition and a one-way hysteresis gap 
of width exactly M. 

To prove this assertion, start the global model mov- 
ing, and let it reach a steady-state at a nonzero velocity. 
Whenever a site hops in the steady state, its local en- 
vironment has necessarily changed on the previous time 
step- cither it has itself hopped on the previous step, 
or one of it's nearest neighbors have. Conversely, any 
site whose nearest neighbor hasn't hopped cannot hop, 
despite the global stress overshoot. So we may apply 
the overshoot stress only to the sites whose environment 
has changed- to those sites that hopped on the previous 
step and their nearest neighbors- with no change in the 
dynamics, and this is the definition of the nonadditive 
model. This completes the proof. 

The nonadditive overshoot model is a local irrelevant 
perturbation of the M — model, which introduces pre- 
cisely M units of one-way hysteresis and does nothing 
else. It does not change the order of the transition or the 
critical exponents. 

The subtraction procedure is inducted from the two 
models defined above. To find a model with no one-way 
hysteresis from a given model, remove the non-additive 
component of the stress overshoots. The nonadditive 
component is the minimum value of stress overshoot on 
any of the sites that are capable of hopping. Since mo- 
tion propagates locally from site to site, any model in 
which the stress overshoots are positive and propagate 
faster than the elasticity will have one-way hysteresis of 



the type described. For our specific model at nonzero M, 
we produce the subtracted model by subtracting the non- 
additive component of the stress overshoot. In the sub- 
tracted model, the overshoot stress 5*(x) is zero if zero 
or one neighbor hops, and S(x) = n — 1 if n neighbors 
of a site hop (recall that a site is its own neighbor). To 
recover our original model, we must add S^x)- a non- 
additive stress overshoot. In other words, S2 (x) is the 
Boolean OR of all 2d+ 1 neighboring velocities. Because 
in steady-state every site that hops has a neighbor that 
hops, the only effect of S2 is to add one-way hysteresis of 
gap size exactly M [l8| . 

The subtracted model is a local perturbation of the 
M = model, and the perturbation is necessarily irrel- 
evant in a renormalization group sense, since it is local 
and breaks no symmetry. We conclude that any hystere- 
sis in the subtracted model is a finite-size effect in this 
regime, a result confirmed by numerical simulation (see 
Fig. I). 

Any model which propagates motion from one site to 
another and which includes a stess-overshoot will generi- 
cally display one-way hysteresis. The appropriate sub- 
tracted model removes the minimum stress overshoot 
which is applied to any moving site. The remaining per- 
turbation is irrelevant, and does not change the critical 
behavior for small values of M. Since the phenomenon is 
so general, we believe that it explains the origin of one- 
way hysteresis in at least some of the systems that display 
it. This implies that one-way hysteretic transitions are 
still continuous, despite appearances. 

Simulating our original model and the subtracted 
model at large values of the parameter M, we find 
that there is a tricritical point at a finite value M c of 
M, beyond which the hysteresis changes character. For 
M > M C) subtraction no longer removes the hysteresis, 
and the transition is truly discontinuous. At M = M c , 
the universality class of the transition is changed. The 
tricritical point is a generic feature of overshoot depin- 
ning models. 

There are two peculiar properties of the hysteresis in 
the unsubtracted models. First, the nucleation radius is 
not infinite as it is in first-order thermal transitions. If 
one well-placed site were to hop when the applied force 
is in the hysteresis gap, the whole surface would start to 
move. Second, the hysteresis loop does not vanish for 
large systems as the inverse logarithm of the system size. 
When the surface is static, the probability of nucleation 
is zero because there is no dynamics. If we were to add 
a small amount of spontaneous hopping, the hysteresis 
loop would vanish as a power of the system size and the 
time. In first-order transitions, the nucleation bubbles 
becomes larger and exponentially more improbable as the 
critical point is approached; here, a nucleation bubble 
may be microscopic. It stays microscopic even at the 
critical transition. 

The characteristics of one-way hysteresis are: 
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FIG. 1: Plot of the finite-size hysteresis gap for the stress 
overshoot and subtracted models. To generate the data, we 
start the system in the moving phase and lower F until the 
surface stops. We then increase the F back up again until the 
surface begins to move. The symbols denote the data, and the 
curves denote the least-squares fits. The error-bars indicate 
la statistical uncertainty. When M = 0.0, we expect the gap 
to vanish as a power law in the system length; the exponent 
is 1/v — 4/3 to first-order in perturbation theory [13,14], 
where v is the correlation length exponent. We measure 1/v — 
1.29 ± 0.03. For the subtracted model, when M = 0.4, 1/v = 
1.26 ± 0.03 and for M = 0.8, 1/v = 1.27 ± 0.03. With the 
subtracted model, there is an M-dependent coefficient that 
increases with increasing M. For the stress overshoot model 
with M = 0.4, the gap is 0.4 within one a at least for the 
larger sizes. We have multiplied this data by \ for plotting 
purposes. For all data, d = 2 and P — 4 [19] . 



1. The stationary phase is unstable. 

2. The nucleation radius is finite at the transition point 

3. With local fluctuations, the probability of nucleating 
a bubble vanishes as a power of the system size like a 
finite-size effect, not logarithmically like a first-order 
hysteresis loop. 

4. The stable phase is described by a continuum theory 
at the critical point. 

The one-way hysteresis which we have seen occurs 
is, in some sense, a phony hysteresis. The reason is 
that it is only present because when the interface is 
static, all overshoots are absent and do not contribute 
to the motion. If we were to add a small amount of 
spontaneous hopping, the hysteresis gap would vanish 
as a power of the system size. The only reason we 
may observe this effect is because in many depinning 
transitions thermal fluctuations may be ignored. 



An earlier model |]| included stress-overshoots only on 
nearest neighbors, but did not include overshoots on the 
site itself. Because of this, the nonadditive component of 
the stress overshoots is zero. If a small amount of self- 
stress overshoot is added, the system acquires a small 
amount of one-way hysteresis. Without the self-stress, a 
site which hops twice in a row without being triggered 
by its nearest-neighbors on the second hop does not feel 
the stress overshoot the second time around. The site 
then has a chance to land and be pinned very weakly, 
and a weakly pinned site may serve as nucleation cen- 
ters when the force is increased in the static phase. The 
result is that the hysteresis gap vanishes as the system 
size grows, albeit slowly because the nucleation centers 
are so rare. The interesting observation is that the nucle- 
ation centers are one (well-placed) site wide. This is the 
numerical and experimental signature of phony hystere- 
sis. A scattering of microscopic seeds will destroy it. In 
solid-on-solid friction, for example, a sound wave with a 
fixed frequency and amplitude should destroy any phony 
hysteresis, while it would only reduce, but not eliminate, 
the hysteresis in a true first-order transition such as the 
one that occurs for M > M c . 

For completeness, we note that when S*(x) = w(x), 
that is, when a site only applies a stress overshoot to 
itself, there is no hysteresis of any kind for any value 
of M no matter how large. This is an accident of the 
parametrization of the lattice model. If the overshoot 
stress is decays exponentially over many time steps, for 
sufficiently slow decay rates we approach the tricritical 
point in these models as well. There is no one-way hys- 
teresis in this model of course, because the nonadditive 
stress overshoot is zero, since the overshoot stress is not 
applied to the neighbors. 
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